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Chiral random matrix theory makes very detailed predictions for the spectral cor- 
relations of the QCD Dirac operator, both in the bulk of the spectrum and near 
zero virtuality. These predictions have been successfully tested in lattice QCD 
simulations by several groups. Moreover, the domain of validity of random matrix 
theory has been predicted theoretically and identified in lattice data. In this talk, 
the current numerical evidence is reviewed. 



1 Introduction 

We are interested in learning as much as we can about the spectrum of the Eu- 
clidean QCD Dirac operator,^) = $ + ig^. One of the reasons for this interest 
is the Banks-Casher relation, S = \(tpi/>}\ = np(0)/V, which relates the chiral 
condensate to the spectral density, p(X) — (VJ n 5(A — A„)}, of the Dirac op- 
erator at zero virtuality. Thus, the spontaneous breaking of chiral symmetry, 
a nonperturbative phenomenon with profound consequences for the hadron 
spectrum, is encoded in an accumulation of the small Dirac eigenvalues. 

It has been realized that chiral random matrix theory (RMT) is a suit- 
able tool to compute the distribution and the correlations of the small Dirac 
eigenvalues, see the talk by Verbaarschot in these proceedings and references 
therein. The first numerical (evidence in support of this statement came from 
instanton liquid simulations £1 The purpose of this talk is to verify the RMT- 
predictions by comparing them to lattice data. We shall see that the agree- 
ment between theory and numerical experiment is remarkable. 

Recently, it has been showna that the results obtained in RMT can also 
be derived directly from field theory (partially quenched chiral perturbation 
theory). Conceptually, this is a big step forward, while for actual calculations 
RMT seems to be the simpler alternative. It is comforting to know that the 
results of the two approaches agree with each other and with lattice data. Re- 
lated topics are covered in the talks by Damgaard, Akemann, Papp, Markum, 
Stephanov, and Halasz (in order of appearance) in these proceedings. 
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2 Symmetries of the Dirac operator 

Because of {1^,75} = 0, the Dirac operator falls into one of three symmetry 
classes corresponding to the three chiral ensembles of RMT: chiral orthogonal 
(chOE), unitary (chUE), and symplectic (chSE) ensemble, respectively.0 In 
most cases, the Dirac operator does not have additional symmetries and is 
described by the chUE. Exceptions are as follows (N c = number of colors). 

- continuum, N c = 2, fermions in fundamental representation: chOE 

- continuum, fermions in adjoint representation: chSE 

- lattice, N c = 2, staggered fermions in fundamental representation: chSE 

- lattice, staggered fermions in adjoint representation: chOE 

The overlap Dirac operator on the lattice has the symmetries of the continuum 
operator. There are a few exceptions where the Dirac operator does not have 
the usual chiral symmetries and is, therefore, described by the non-chiral RMT 
ensembles: QCD in three dimensions (UE), and the Wilson Dirac operator on 
the lattice (OE for N c — 2, UE for N c > 3, SE for the adjoint representation). 

3 Correlations in the bulk of the spectrum 

The full Dirac spectrum can be computed numerically using, e.g., a special 
version of the Lanczos algorithm. Which features of the spectrum are de- 
scribed by RMT in the bulk of the spectrum, i.e., away from the edges? The 
global spectral density is certainly not given by the RMT result since it is sen- 
sitive to the details of the dynamics. However, one can separate the average 
spectral density from the spectral fluctuationSpOn the scale of the local mean 
level spacing. This process is called unfoldingfl After unfolding, the spectral 
correlation functions are, up to a certain energy, given by RMT. This limiting 
energy is called the Thouless energy and will be addressed in Sec. |^. 

The short- and long-range correlations of the eigenvalues are measured 
by quantities such as the distribution, P(s), of spacings, s, between adjacent 
eigenvalues and the number variance, Y?{L) = ((n(L) — (n(L))) 2 ), where n{L) 
is the number of eigenvalues in an interval of length L. In the bulk of the 
spectrum, the predictions of the chiral RMT ensembles for these quantities 
are identical to those of the corresponding non-chiral ensemble. They can 
be tested against lattice data. The predictions of RMT were confirmed wilt 
very high- accuracy in the following cases: staggered fermions in SU(2),B. 
in SU(3),H and in compact U(1),Q respectively, Wilson fermions in SU(2)J! 
and the overlap Dirac operator fix SU(2), in SU(3), and in the adjoint rep- 
resentation of SU(2), respectively]!] Furthermore, the lattice-xlata continue to 
agree with the RMT predictions in the deconfinement phase.Q In all cases, the 
agreement was perfect, see P(s) in Fig. [I] for an example. 
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There is a very interesting question concerning staggered fermions: for 
N c = 2, they have the symmetries of the chSE whereas the continuum sym- 
metries are those of the chOE (and the other way around in the adjoint 
representation). While one should eventually see a transition from chSE to 
chOE behavior as the lattice spacing goes to zero, it is unlikely that this can 
be observed on present day lattices. The overlap Dirac operator, on the other 
hand, has the correct symmetries in all cases. 

If a chemical potential is added to the problem, the Dirac eigenvalues are 
scattered in the complex plane. In this case, P(s) can also be constructed, 
see the talk by Markum. 

4 Correlations at the spectrum edge 

Because the nonzero eigenvalues of ift come in pairs ±A„, the Dirac spectrum 
has a "hard edge" at A = 0. As mentioned in the introduction, the Dirac 
eigenvalues in this region are of great interest because of their connection 
to chiral symmetry breaking. The distribution of the smallest eigenvalues is 
encoded in the microscopic spectral density, p s (z) = limy^oo p(z/VT,)/VT,, 
which is universal and can be computed in RMT. One can also construct the 
distribution of the smallest eigenvalue, P(A m ; n ), and higher-order spectral cor- 
relation functions. Analytical RMT results for these quantities are available. 
To compare these to lattice data, the only input one needs is the energy scale 
VY, = Tvp(0), which is obtained from the data by extracting p(0). 

The spectral correlations in the vicinity of A = 0, in contrast to those in 
the bulk, are sensitive to the number of massless (or very light) quarks and 
to the topological charge v. This is another reason why the hard edge of the 
spectrum is more interesting than the bulk. 

4-1 Quenched approximation 

The first evidence from the lattice for the universality of the microscopic 
spectral density came from an analysis by Verbaarschot of Columbia data 
for the valence quark mass dependence of the chiral condensate in SU(3) 
with staggered fermions.B Plotting all the data computed for various values of 
/3 = 6/g 2 in a rescaled form suggested by RMT, it was observed that in the 
region of very small masses they all fell on the same universal curve given by 
RMT. The simulations were done with dynamical quarks, but the sea quarks 
were much heavier than the valence quark so that one effectively had Nf = 0. 

The first direct lattice computation of the microscopic spectral quantities 
was done in quenched SU(2) with staggered fermionsjlj and the data were 
subsequently analyzed in more detaiLEU The agreement between lattice data 
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Figure 1. Spectral properties of the QCD lattice Dirac operator. The histograms represent 
lattice data, the curves are RMT predictions. Left: Nearest neighbor spacing distribution 
for SU(3), staggered fermions (chUE), V = 6 3 X 4, (3 = 5.0H Middle: Microscopic spectral 
density for SU(2), staggered fermions (chSE), V = 10 4 , = 2.oJl3 Right: Distribution erf 
the smallest eigenvalue for SU(2), overlap Dirac operator (chOE), V = 4 4 , j3 = 1.8, v = 1.0 



and RMT predictions was remarkably good, see p s (z) in Fig. [I] for an example, 
but the analysis immediately raised a number of questions. First, why were 
the data consistent with the RMT predictions for topological charge v = 
even in weak coupling? Second, at what energies does RMT cease to be 
applicable? The answers will be given in Sees. |J and §. 

Meanwhile, several groups have performed lattice simulations of the hard 
edge of the spectmm for almost all interesting cases: SU(3) with staggered 
fermions in three&3 and fouiil3llil dimensions, theJ3chwinger model using both 
a fixed point and the Neuberger Dirac operator Jij SU(2) and SU(3) with stag- 
gered fermions in the adjoint representation,!^ and the overlap Dirac operator, 
in SU(2), in SU(3), and in the adjoint representation of SU(2), respectively!] 
The theoretical predictions of all three chiral ensembles (and of the UE) have 
thus been verified with high accuracy. 



4-2 Light dynamical fermions 

In the chiral limit, the microscopic spectral correlations predicted by RMT 
depend on the number, Nt, of massless dynamical quarks. If the dynamical 
quarks are sufficiently heavy, the microscopic spectral correlations are given 
by the quenched result. There is an intermediate regime for quark masses of 
order 1/VE, i.e., of the same order as the low-lying Dirac eigenvalues. In this 
"double-scaling" regime, the RMT results depend on the quark masses. So 
far, analytical results are only known for the chUE. Lattice simulations with 
very light-dynamical quarks have been performed in SU(2) with staggered 
fermions, Ell corresponding to the chSE for which the corresponding RMT pre- 
dictions were computed numerically. Again, very good agreement was found 
for p s (z) and P(A mi „). 
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4-3 Topology 

Lattice simulations with staggered fermions show that the microscopic spec- 
tral quantities agree with the RMT predictions in the sector of zero topological 
charge, even in weak coupling. Presumably, the reason is that the would-be 
zero modes related to topology are shifted away from zero due to discretiza- 
tion errors of order a 2 , where a is the lattice spacing. Thus, it is essential to 
use "better" Dirac operators which obey the Ginsparg- Wilson relation and for 
which the Atiyah-Singer index theorem holds. This has recently been done, 
in the Sciwvinger model using both a fixed point and the Neuberger Dirac 
operator,t3 and using the overlap Dirac operator in SU(2), in SU(3), and in 
the adjoint representation of SU(2), respectively.^ ("Neuberger" and "over- 
lap" refer to the same operator.) The lattice data agree very well with the 
RMT predictions also in sectors with v ^ 0, see P(X min ) in Fig. [I]. 

5 Thouless energy 

Since RMT predictions are derived without any knowledge of the details of the 
QCD dynamics, they are only valid up to a certain energy which is called the 
Thouless energy, E c . To identify E c for QCD, one needs two ingredients: (1) 
the fact that the random matrix model is applicable only if the kinetic terms 
in the chiral Lagrangian can be neglected, which is the case for L <C l/m r 
(with L the linear extent of the four- volume and the pion mass), and 
(2) the Gell-Mann-Oakes-Renner relation, connecting to the quark mass 
and the chiral condensate. This gives rise to the theoretical predictior£3o 
E c ~ / 2 /I]L 2 , where f„ is the pion decay constant. This prediction was 
verified numerically in inptanton liquid simulationalS and in lattice simulations 
for SU(2)E2rQ and SU(3)E-3 with staggered fermions. For lattice simulations, 
this roughly means that the domain of validity of the RMT approach is larger 
for larger lattice size and stronger coupling. 

6 Summary 

We now have a wealth of analytical and numerical evidence in support of the 
claim that the low-lying spectrum of the QCD Dirac operator is described 
by universal functions which can be computed, e.g., in RMT. This applies to 
the quenched approximation, to the full theory with dynamical quarks, and to 
trivial and nontrivial topological sectors. Moreover, we know the energy where 
the universality breaks down. Apart from a better analytical understanding of 
the Dirac spectrum, this also offers a variety of practical applications, such as 
better extrapolations to the thermodynamic limit £i the extraction of chiral 
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logarithms,E3 and perhaps the construction of hybrid Monte Carlo algorithms 
which would use the distribution of the small eigenvalues as input. 
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